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1. Introduction 

During the last decades the theoretical investigation of strongly correlated electron 
systems concentrated on the one-band Hubbard model |], which serves as the standard 
model for electrons with strong short-range interactions; for a recent review, see |2j . The 
Hubbard model was originally proposed for the description of 3d electrons in transition 
metals to explain ferromagnetism and antiferromagnetism in iron and nickel and their 
oxides |1|, |3|]. For these substances the band degeneracy and intra-atomic (Hund's rule) 
exchange couplings apparently play an important role and, consequently, multi-band 
Hubbard models need to be investigated. 

Until recently, band degeneracies were considered mostly for immobile / electrons 
which hybridize with featureless conduction electrons. The corresponding single- 
impurity Anderson model f|] and its periodic generalization (Varma-Yafet model ||) 
apply, e.g., for the rare-earth materials. Since the Coulomb interaction between the 
/ electrons is very large, one often encounters the situation where the iVy-fold degenerate 
levels are at most singly occupied. Hence, the issue of degenerate bands with finite 
electron density need not be addressed in these cases. 

In the past few years the multi-band Hubbard model received new attention. 
The Gutzwiller approximation to Gutzwiller-correlated variational wave functions was 
generalized in @-[f§, the Slave Boson mean-field approach along the lines of Kotliar 
and Ruckenstein was developed and used in flOfl-ll^l, and dynamical mean-field 



methods were applied in [TJ|, 14]. In this paper we extend the method in [15, ljj 



to evaluate general Gutzwiller-correlated wave functions in the limit of infinite space 
dimensions without further approximations. In this way we recover the results of |J-[H] 
and thus show that these earlier results become exact in infinite dimensions within the 
variational approach. Furthermore, we prove that the Slave Boson mean-field results 
are variationally controlled at zero temperature in this limit. Similar results hold for 
the one-band Hubbard model; see |T7J and 0, Chap. 3.4, 3.5] for brief reviews. 



Our paper is organized as follows. In Section g we introduce a general class of 
Gutzwiller-correlated wave functions which allows for the variational study of Hubbard 
models with general spin-orbit structure. We set up a diagrammatic perturbation theory 
to calculate expectation values for the variational ground-state energy and other physical 
quantities of interest. In Section || we use our formalism to derive exact analytical 
expressions for the ground-state energy in the limit of infinite spatial dimensions. In 
Section |4] we compare them to the results of approximate treatments of Gutzwiller- 
correlated wave functions and those of the Slave Boson mean-field theory. Since our 
treatment allows for the inclusion of site-diagonal energetic disorder in the Hamiltonian, 
we briefly discuss the consequences of strong correlations on the Anderson transition 
within our variational description. A summary in Section |5| closes our presentation. 
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2. Gutz wilier-correlated wave functions 

In this section we introduce the Hubbard Hamiltonian for degenerate bands and the 
class of Gutzwiller-correlated wave functions as our approximate ground states. The 
minimization of the ground-state energy fixes the variational parameters contained in 
these trial states. For the evaluation of the corresponding many-body expectation values 
we develop a diagrammatic formalism which allows for the complete solution of the 
problem in infinite space dimensions. 

2.1. Definitions 

In his first two papers on narrow-band electron systems Hubbard [I], |3| considered 
3c? electrons with a purely local interaction. To simplify our considerations we make 
the further assumption that the interaction depends only on the number densities of 
electrons in the d orbitals. Then we may write the Hamiltonian in the form 

where a, a' = 1, . . . , 2N are combined spin-orbit indices (N = 5 for d electrons), i, j 
denote lattice sites, and t^J is the matrix element for the electron transfer between two 
sites i and j with spin-orbit indices a and a', respectively. In the following we will use 
the notion "orbital" for spin-orbital states. 

The local energies t°f = S a ye^ can be arbitrarily chosen to describe different orbital 
energy levels, e.g., in compounds or to mimic the influence of impurities and other 
lattice defects. As we will show below, expectation values for Gutzwiller-correlated 
wave functions can be calculated analytically in infinite dimensions even in the presence 
of energetically random impurity potentials. 

Finally, U i ' describes the local Coulomb interaction between two electrons in the 
orbitals a and a' on the same lattice site. The local interaction partly accounts for the 
atomic exchange coupling. For example, in the case N = 2 the ground state of an atom 
with two electrons should be a spin triplet according to Hund's first rule. According 
to our Hamiltonian ([]]) the interaction energy of the local S z = ±1 states is different 
from that of the S z = states. However, the two S z = states are still energetically 
degenerate in (|l|), i.e., the interaction does not fully distinguish between triplet and 
singlet states, as required by Hund's first rule. In this work we restrict ourselves 
to density-dependent terms only. This should be the generic case for C 60 and other 
fullerenes for which exchange coupling is less important because of the large diameter 
of the molecules; see, e.g., |18| . For other materials it might be necessary to include 
local spin-flip terms both in the Hamiltonian and in the variational description. Since 
the evaluation of many-body wave functions is fairly complicated we restrict ourselves 
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to Gutzwiller correlators which are solely dependent on the density operators [see (7) 
below], and also omit the spin- flip terms in the Hamiltonian for consistency. 

We introduce the following notations for the 2 2N possible configurations of a given 
lattice site. 

(i) An atomic configuration / is characterized by the electron occupation of the 
orbitals, 

J6{0;(l),...,(2iV);(l,2),...,(2,3),...(2iV-l,2iV);...;(l,...,2iV)} , (2) 

where the order of numbers in (abc. . .) is irrelevant. The symbol in (0) means 
that the site is empty. In general, we interpret the indices I in (||) as sets in the 
usual sense. For example, in the atomic configuration I\I' only those orbitals in / 
are occupied which are not in I'. The complement of / is I = (1, 2, ... , 2N)\I, i.e., 
in the atomic configuration / all orbitals but those in / are occupied. 

(ii) The absolute value |/| of a configuration is the number of elements in it, i.e., 

|0| = 0; | (a) | = 1; \(a, b)\ = 2; . . . ; . . . , 2N)\ = 2N . (3) 

(iii) The operator which projects onto a specific configuration I on site i is given by 

rrii-i = Y[ n i;(7 J[ (1 - n i;o ) , m i;/ = (m i;I ) , (4a) 

where (...) denotes the expectation value in the Gutzwiller-correlated wave 
function; see below. The operators m i; / measure the "net" occupancy. Besides 
these we define the operators for the "gross" occupancy as 

n i; i = I| hi;* ! n i;i = (*V> ! ^;0 = 1 ■ ( 4& ) 

crel 

The gross occupancy operator h^i gives a non-zero result when applied to /' only 
if / contains electrons in the same orbitals as However, / and /' need not be 
identical because I' could contain additional electrons in further orbitals, i.e., only 
/ C I' is required. 

Each gross (net) operator can be written as a sum of (net) gross operators 

Hi = H ' ( 5a ) 

I'DI 

mj= £(-l) |/VI "i;/'- (56) 

I'DI 

For practical calculations the net operators vn^i are more useful than the gross 
operators fii-j because the former are projection operators onto a given configuration i" 
on site i, i.e., rh^i rhi-ji = Sij'Thi-j. 
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With these definitions we may rewrite the interaction part of the Hamiltonian ([J) 



as 



E U^ a hi-vhi-v' = E Ui-jrhi-j , (6a) 

i;a,a' i;I (\I\>2) 



with 



Ui-j = E u r' • (66) 

a,a'el 

With the help of these definitions we may formulate the class of Gutzwiller-correlated 
wave functions for degenerate bands as 

|*g> = Pq |*o> , (7a) 

Pg = n n £r ■ m 

i I 

(m>2) 

The trial states depend on = 2 2N — (2N+1) real numbers g^i for each lattice site and 
on further variational parameters in \^q). I n general, the gross occupation densities n i;o - 
in |^g) are different from those in \^q). The generalized Gutzwiller correlator Pq in (7) 
suppresses fluctuations in the multiple orbital occupancy for repulsive interactions. In 
a translationally invariant system the Gutzwiller wave functions for degenerate bands 
as proposed in |§ are recovered by setting g^i = gj. Similar but different expressions 
for the Gutzwiller wave function for degenerate bands can be found in the work by 



Gutzwiller [IS], Gutzwiller and Chao [20], and Chao [21 



2.2. Diagrammatic evaluation 

To gain further insight into the physics of the variational wave functions we have to 
evaluate expectation values 



(6) = \ri r . (8) 





0|*G> 


(*g| 





The variational parameters in |^g) are obtained by the minimization of the expectation 
value of the Hamiltonian (|l|), 

E™= min(F), (9a) 
(H) = tifictA-y) + E U,j(m,j). (96) 

i,j;cr,a> i;I (]J]>2) 

The variational ground-state energy E™ is an upper bound for the exact ground-state 
energy This upper-bound property applies only if we are able to evaluate the variational 
ground-state energy without further approximations. 
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The evaluation of expectation values with correlated wave functions is a many- 
particle problem that cannot be solved in general; see |17]] and || Chap. 3.4] for a review. 



In this paper we will use the method introduced by Gebhard |15|, |16| , which allows the 
approximation-free evaluation of general Gutzwiller-correlated wave functions in the 
limit of infinite spatial dimensions. The formalism is based on a diagrammatic expansion 
of expectation values in such a way that the (variational) self-energy identically vanishes 
in infinite dimensions. 

As shown in more detail in [[OJ we have to carry out the following program: (i) choose 



the appropriate expansion parameters, (ii) apply Wick's theorem, and (iii) use the 
linked-cluster theorem. If the expansion parameters are properly chosen the lowest 
order in the expansion already gives the exact result in infinite dimensions. 

2.2.1. Choice of the expansion parameters. As a first step we choose the appropriate 
expansion parameter(s) for our diagrammatic theory. In this respect the variational 
approach is more flexible than the standard perturbation theory for Green functions in 
interacting electron systems. 

We express the one-particle product wave function \^o) in the form 

i^o)=riwi!i i ---^i < fo) , (io) 

i 

where |$o) is another, normalized one-particle product wave function. The real numbers 
Vt$i ■ ■ ■ iVipN are chosen such that the square of the "modified" (local) Gutzwiller 
correlator, 

27V 

(7=1 I 

tm>2) 

2V 



= T7«;0 II [1 + (%r _ II [1 + (9i;I ~ , (H) 

a=l 

can be written in the form 



CT=1 I 

tm>2) 



Pi = 1 + 2^1,2^2 + ^1,37^3 + • • • + ^;l,...,27V^i F ... i 2JV ■ ( i2 ) 

Here, the parameters x^j are real numbers and the Hartree-Fock (HF) operators h^J 
are defined as 

^?=II(V-0 ^ |/|>1, (13) 

and 7T,?J = 1. Note that the definition of the Gutzwiller correlator and the Hartree-Fock 
operators differs from the one given in [15, [16| . In (|13|) we introduced 



($ |^;cr|$o) = (n i; a)o (14) 
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as the local densities in orbital a in the new single-particle product state |$o)- 
Equation (|i~2^) poses 2 2N conditions for the 2 2N — (2N + 1) parameters x^i (|/| > 2) and 
the 2N + 1 parameters 77^0, . . . , t^tv- We will solve for them in terms of the original 
variational parameters g^j > 2) in the next section. 

The parameters x^i (|/| > 2) go to zero for small interaction strengths since 
9i-,i(Ui ,cr — ► 0) — ► 1. Hence, we may use them as the expansion parameters for 
a perturbative approach. In its diagrammatic formulation the Xi-i play the role of 
(internal) vertices at which (at least) four lines intersect. The crucial point in the 
expansion flT2"|) is the fact that there will be no Hartree ( "bubble" ) diagrams which are 
of order unity in all dimensions. Since we are interested in simple expressions in infinite 
dimensions, we included their contribution in the expansion parameters Xi-j. Note that 
we will make the necessary assumption that local "Fock terms" do not occur, i.e., we 
demand 

($o|c+ CT c i;CT ,|$ ) = <V n i;a (15) 
for the one-particle product wave wave function |$o)- 

2.2.2. Application of Wick's theorem. As our second step we formally expand the 
expectation values that we need for the calculation of the variational ground-state 



energy (9). With the help of (12) we may write 

(* G |m /;/ |* G ) = ($ \P f m f .jP f J] P-\%) , (16a) 

¥/ 

(* G |C+ CT/ C M J* G > = ($ | (PfSf^Pf) (PhC h ;aA) U P "\^), (166) 

&f(f,h) 

n^ = i+£^E' e n (^41) . (i6c) 

i^f,h k=l ' %i,—,ik I i 1 >—> I i k j=ii 

i\h\>2) 

where the primes on a sum indicate that all lattice sites are different, 

f±h±h±...±i k . (16d) 

There are still some Hartree contributions contained in (16 a) and (166) which come from 
the "external" sites / and h. With the help of (4) we can always find unique expansions 
of the form 

if 

PfcU P f = d U E (176) 

If{aj0f) 

hc h , a P h = c h .^ £ zZAl- ( 17c ) 



Later, in Sect. |3^, we will calculate explicitly the coefficients oj. If and Zj! If in terms of 
the variational parameters g^i and 77^0, T]i ;a . We introduce the abbreviations 



n f;l f 11 "./:/, 
j=h 



0/ > 



4£VV/. <*) = (*«> 



C f;a f C h;a h n f 



iHF r? HF TT r? HF 

;// ra ft;/h 11 n j;/. 



0/ > 



'18a) 



;i86) 



3=n 



where the products are replaced by unity for k — 0. This notation allows us to 
rewrite (16) as 



(* G |m /; /l*G> = E°k 



fe=l ' ii,...,i k rj 1 ,.-.,li fc \j'=n 



,1k 



(19a) 



X 



CO 1 2fc 

^V/^o+Ei? E' e (n^0^i:* /lfc K,< 

fc=l ' »i,...,ifc h 1 '---'h k j=%i 

(l^l>2) 



(196) 



The primes on the lattice sums indicate that all lattice indices are different from each 
other. Moreover, the definition (||) implies that a given orbital index occurs only once 
at each lattice site. 



Now we are in the position to apply Wick's theorem [^] . The expectation values (18) 
can then be written as determinants. Equation (18a) becomes 



i f ,i ilt ...,ii k 

f,h,—,ik 



Mf, f M fM ■■■ M Uh 
M hJ M luil ••• Mi 



M ih j M ihtil ■■■ Mi 



ik,ik 



(20) 



where we introduced the sub-matrices 



Mi 



1,3 



'■J 



I p a \k\> a i 



p >'3> 



(21) 



/ 



Here, the orbital indices <7i, . . . , a\i t \ (a[, . . . ,aL<) are the elements of U (Ij) and 



p°f = (i - ^x^ol^a^i^o) = (i - kMA^ ( 22 ) 

is the one-particle density matrix for the one-particle product wave function |$o) for 
% ^ j. These objects play the role of "lines" in our diagrammatic expansion. Note that 
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we do not have to distinguish between "hole" and "particle" lines because all sites are 



different when we apply Wick's theorem [17|. Furthermore, local terms did not arise 
because we subtracted the Hartree contributions and ruled out Fock terms according 
to (|15"D. Similarly, equation (186) becomes 



S 



f,h,ii,...,ik 







M h>h 


M hM ■ 


■■ M h,i k 






Mf,K 


Mf, h ■ 


■■ M h,i k 


_l\\If\\Ih\-\If\-\Ih\ 




M h)h 


M luil • 


■■ M hjik 




M ikJ 






■■ M ik>ik 



(23) 



The matrices M it j in ( p3|) are again given by (|2l|) . Here, the orbital indices belonging 
to the lattice sites / and h are elements of //, // U 07, I h , and Ih U cr h , respectively. 
To calculate expectation values for \^g) we have to divide the numerators, (19), by 



the norm 



OO 1 Ik 

(*g|*g> = <*oi*o> + e F e' e (n ^) T t:t 

k=l ' h,—,ik I i 1 ----h k j=n 

(|iil>2) 



(24) 



In principle, we could derive diagram rules for the series expansion of the determinants 
in powers of the parameters Xj.j.. In this paper we restrict ourselves to the limit of 
infinite dimensions where not a single diagram needs to be calculated. 



2.2.3. Linked- cluster theorem. The summation restrictions prevent us from the 
application of the linked-cluster theorem. In the case of a single band (N = 1) the 
summation restrictions (16 d) can simply be dropped because the determinants (p0|) 
and fl2"3"|) vanish identically if two lattice indices coincide. Since we may then 



independently sum over all lattice sites the linked-cluster theorem applies such 
that the disconnected diagrams in the numerator are canceled by the norm ||15|| . 

The case N > 2 requires more care because, in general, the determinants T in (p0|) 
and S in (|23D remain finite when we equate two lattice sites. To make progress we note 
that a summation restriction over the spin-orbit indices in I ix , . . . , I ik will not prevent 
the applicability of the linked-cluster theorem. To see this, we write a typical spin-orbit 



sum in (19) in the form |L6| 

2N i 

E/W = Ei7 E' /(«!>•• 



, av 



r=l ' " oil, 



(25) 



The ensuing summation restriction «i ^ «2 ^ • • • ^ « r can be removed because the 
determinants in ( p0|) and (|23|) vanish if two of these indices are identified. 

Now we equate two lattice sites / and m, e.g., in the determinant (|20|). Two different 
cases have to be distinguished: (i) if J/ and I m have at least two elements in common 
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the determinant (|20|) vanishes since at least two rows or columns will be identical; 
(ii) otherwise, if J/ and I m have no common element, the determinant (|20|) remains 
finite. When we equate two lattice sites in the determinant we effectively map two 
vertices onto each other in the corresponding diagrammatic expansion, i.e., we generate 
a diagram that already appeared at some lower order in the expansion; see [16[ for 



a simple example for this "vertex packing". The original diagram with I and m put 
equal and the corresponding lower-order diagram have the same topology but different 
prefactors, Xi-j lU j m for the lower-order diagram and Xi-j l x rn .j m for the new contribution. If 
we relax all summation restrictions in the numerator and in the norm we must replace 
the vertices x i; i. and the coefficients of. 7 . and zf.j by effective ones, i.e., x^j — ► x i; i, 
~ > an d z i-,i ~^ z ti- Fortunately, the lowest-order coefficients o\.$ and z£.q 

remain unchanged by this procedure because these coefficients correspond to diagrams 
without an internal vertex whereas the "vertex packing" leaves behind at least one 
internal vertex fT6j. In the next section we will show that o\.§ and zf^ alone give the 



exact result in infinite dimensions. Since the number of correction terms generated in a 
given order of the expansion remains finite, the calculation of systematic 1/d corrections 
is still possible although very tedious. 

Now that we formally eliminated all summation restrictions we can apply the linked- 
cluster theorem. We thus find 

oo I r *fe ~| C 

(rhf.j) = 0% + J2 5l f;i f E u E E Kf-Jf II -0:/. ".!!/'. f ( 26 «) 

If k=l ^ h,...,i k l iv -,I ih L j=h J 

(\If\>v (\h\>2) 

(Cf;a f Ch;a h ) = Z f;$ Z h^l(^f^f^h-,cr h )o + E Z f-If Z h;I h 

(\I f \+\I h \>V 

k=l h,...,i k h v -,h k k j=ii J 

(|Iil>2) 

where, as usual, {. . .}q indicates that only the connected diagrams need to be considered. 



3. Exact results in infinite dimensions 

In this section we briefly review the diagrammatic simplifications which occur in 
infinite dimensions. We find that in our approach not a single diagram needs to be 
calculated in this limit. Consequently, we derive explicit analytic expressions for the 
variational ground-state energy, the one-particle density matrices, and the average net 
occupation densities which are exact in infinite dimensions and valid for the whole class 
of Gutzwiller-correlated wave functions. 
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3.1. Simplifications 

Systematic studies of the limit of infinite dimensions for itinerant electron systems 
started with the work of Metzner and Vollhardt p5[ ; for details and a recent review, 
see [fj, Chap. 5]. One of the essential simplifications is the following: if two vertices I 
and m are connected by three independent (Green function) lines only the contribution 
for I = m survives in infinite dimensions. For example, in the one-band case the 
(proper) self-energy becomes purely local in this limit. In our variational theory the 
lines between two sites represent the one-particle density matrices Pf£ (|22|) . They 
vanish by construction, if we set I = m. Consequently, the variational (proper) self- 
energy is identically zero in the limit of infinite dimensions in the one-band case [15|. [T(| . 
Note that we guaranteed in our expansion that there are no Hartree contributions such 
that all (internal) vertices are connected by at least three lines. For the case N > 2 
vertices with more than four lines appear for which our arguments particularly apply, 
and any diagram with more than one line must vanish in infinite dimensions. 

Consequently, not a single diagram needs to be calculated. Instead, we immediately 
find from (26) 

(rhf.j)^ =0%, (27 a) 

in infinite dimensions. Thus, the problem is formally solved since only the coefficients 
0^.0 and Zf f Q and the properties of the single-particle product wave function |$ ) enter 
the final expressions. 



3.2. Explicit results for the ground-state energy 

Thus far the variational ground-state energy for general Gutzwiller-correlated wave 



functions could only be accomplished for N — 1 [I5[ [T7[| . Here we show that explicit 
expressions for the ground-state energy can be obtained for all interaction strengths and 
all AT > 1. In the following considerations we suppress the spatial index. 

The remaining problem is the calculation of the coefficients o\.§ and z£.q. We recast 
the modified Gutzwiller-correlator (11) into the form 



m 



(2f 



where we introduced the notation 
fj$ = 1 , 

Vi = 
9i = 



II ^ for |/| > 1 
1 for \I\ < 1 



(29) 



12 

Then, the left-hand side of (17a) becomes 
PrrtiP = rjlvfal™! ■ 



(30) 



Now we expand the net occupancy operators mj in terms of the Hartree-Fock operators 



* HF 

nf as 



mi 



— n 



HP 



zk-rr 7 ' 1 n n° n 0-<) 



(31) 



We compare this expression with the right-hand side of (17a) and find the simple result 



4 = (-i) |/nj,| %w n ^ n 

<tG/\/' <re/\/' 

The representation of the modified Gutzwiller correlator in (^) allows us to write 

L / 

Again, we expand the net occupancy operators rhi in terms of the Hartree-Fock 
operators fif F (31). 

Now we are in the position to compare ([33]) with (O). The constant coefficient gives 



(32) 



(33) 



where 



m 



1 + E {vWi - l) m] 



?=n«°n(i-!) 



(34a) 



(346) 



is the corresponding net occupation density in the uncorrelated wave function |$o)- By 
construction, the coefficient to first order in n^ F is zero in (1121). We thus find 



(35) 
fix the 



nl E E («M - 1 r^o = ° ■ 

Finally, the higher orders in the Hartree-Fock operators in (^) and 
coefficients £/, 



~2 2 ' 



(36) 



where we used (|32]) . For N > 2 one cannot explicitly solve (34), (pop, and ( |36| ) for 
r/o-, and X/ in terms of n° and gj. 
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In the limit of infinite dimensions we can obtain explicit results for the ground-state 
energy. From (27a) we find 

mi = (mi) = o\ = T]^g 2 I m° I (37) 

in infinite dimensions. Equation fl3"7|) implies that the ?7-terms can be expressed by the 
net occupancy densities as 

2 m 2 m a m% 

Note that these simple results do not hold in terms of the single-particle product wave 
function |\& ) in (7a) but only in terms of |$ ); compare equation fllpp. Furthermore, 



equation p7| ) allows us to replace the original variational parameters gi by their physical 
counterparts, the net occupancy densities mj > 2), as 
2 _ m ab m$ 
9a > b ~ m a m b ' 

: (39) 

9\,...,2N ~ 

mi---m 2 N 

Equations (39) are well known from the Gutzwiller approximation for the one- 
band |25], |3 an( i the multi-band case |§. They show that the parameters gj rule 



the law-of-mass action between single occupancies of a site on the one hand and its 
multiple occupancies and vacancies on the other hand. 

In infinite dimensions equation (34) is trivially fulfilled since 

m0 = 1 — mj (40) 

is true by definition. From the definition of net and gross operators (4) we know that 
the (gross) occupancy in the orbital a is given by 

n a = m i ■ ( 41a ) 
We further use 

^ 771/ = 1 — 7710 — n a (41 b) 

7|>1 



in (|35|) . With the help of ([38|) it is easy to show that 

n a = 7i° (42) 

holds in infinite dimensions. It is thus seen that the local densities in the Gutzwiller- 
correlated wave function |^g) an d in the one-particle product wave function |$ ) are 
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the same. This had not been so if we had worked with the one-particle product wave 
function |^o)j see (7a). For this reason we introduced the new one-particle state |$o) 
in (|T0|), where the single-particle r\ terms could be interpreted as "chemical potentials" 
which guarantee that the average single-orbital occupancies remain unchanged by the 
modified Gutzwiller correlator (11). 

Thus far we replaced the original variational parameters gi by their physical 
counterparts, the average multiple-occupancies m/. As a last step we have to express 
the coefficients z* in terms of our new variational parameters. According to (96) 
and (27b) they can be interpreted as site-dependent renormalization factors for the 
electron transfer between two sites. For their derivation we introduce the operators 

mf= I] *V I] (1 - n<r>) ■ (43) 

<r'eI\<T cr'el\cr 

Then we may write 

Pc+P = r]Al+ E (Vi9i - ) ctr]Al + E fej-l)™/ 

= vlctU+ E (Wiu ( -l)A?)(l+ E (TO-l)m?) (44) 

= nlK ( 1 + E (viuvgiuaVigi - 1) mj 

With the help of (31) we may again expand the projection operators in terms of the 
Hartree-Fock operators fif F . A comparison with the definition (176) then gives 

Zi=V&i,$ + vl E (-lf nTl (Vi'uvgi>u*r)i'9r-l) U <> II (l - ^) • (45) 

This includes the special case of / = which we need for the calculation of the 
renormalization factors for the electron transfers, 

4 = r ll + vl E (Viu*giuvrji9i - 1) : _ J • (46) 

We use fl37|) and the identities 

71° 

E m? = l-n°, (476) 

I(cr£/) 

see (346) and (41a), and finally arrive at 

V^=4= f, 1 n 51 V^W^f • (48) 
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This leads to our final expression for the ground-state energy 

i,j;<r,a' i;I(\I\>2) 

To fix the variational parameters my (|/| > 2) and those in |$o) the minimum of 
this expression has to be determined. Applications of the formalism are planned to be 
published elsewhere. 



4. Comparison with other approaches 

In this section we compare our results to those of related work, namely the Gutzwiller 
approximation and the Slave Boson mean-field approach, which are seen to become 
variational controlled in infinite dimensions. Finally, we briefly discuss the applicability 
of the variational approach to the Anderson transition in strongly correlated electron 
systems. 



4-1. Generalized Gutzwiller Approximations 

For the one-band case it is known that the results of the Gutzwiller approximation 



for the Gutzwiller wave function become exact in infinite dimensions 15, [23l 26 



In addition, the variational calculation of correlation functions is straightforward in 



infinite dimensions [15], whereas their consistent treatment is beyond the Gutzwiller 
approximation scheme. 

Generalizations of the Gutzwiller approximation to other situations than the 
translationally invariant, paramagnetic one-band Gutzwiller wave function encounter 
two conceptual difficulties, (i) It is not clear from the beginning whether the generalized 
Gutzwiller approximations in lead to physically acceptable results in the whole 

parameter space. For example, early extensions of the Gutzwiller approximation to 



the case of antiferromagnetism gave rise to negative occupation densities ||27|| . (ii) The 
Gutzwiller approximation may give the correct results in infinite dimensions but the 
corresponding wave function could not be identified properly. For example, in a recent 
treatment Okabe found the correct expression for the ground-state energy ( |49l) in 
infinite dimensions for the translationally invariant case but he did not specify the 
corresponding Gutzwiller wave function. Therefore, his results could not be tested 
against numerical approaches. In contrast, the wave function (7) provides a solid starting 
point for variational Monte-Carlo simulations. Hence, an assessment of the quality of the 
Gutzwiller approximation for finite dimensions can now be performed; see, e.g., |28|, ^ 
for applications to the one-band case. 

The proper choice of the wave function (7) was essential for its approximation- 
free evaluation in infinite dimensions. In our definition of generalized Gutzwiller- 
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correlated wave functions each multiple occupancy of a lattice site is controlled by 
its own variational parameter. Equations (39) show that we may then replace these 
variational parameters by their physical counterparts, the net multiple occupancies mj. 
In contrast, the trial state of Chao and Gutzwiller |20| and Chao |21[] only includes 



variational parameters for the gross occupancies nj for |/| = 2, 

i*G>=n n ^ ;/ i*o>, (so) 

i I 

cm=2) 

where we used our notation. In their evaluation the "maximum term conditions" of the 
Gutzwiller approximation also lead to a set of 2 2N — (2N + 1) equations similar 



to (39) but it cannot be solved explicitly because there are only N(2N — 1) variational 
parameters gj (|/| = 2) in fl50|) . To make progress Chao set to zero all terms with 
multiple occupancies |/| > 3, an assumption not warranted by the wave function (|50|). 
If we (artificially) put m i: / = for |J| > 3 for our wave function (7) we recover Chao's 
results. Hence, Chao's results do apply to the Gutzwiller-correlated wave function (7) 
in a special limit. 

Similarly, the results of Lu 0] are found to be correct if his further assumptions are 
adopted. He used site and orbital-independent Coulomb interactions, U[' a = U, and 
assumed that some of the multiple occupancies vanish identically. In [[| it was shown 
that the latter assumption might not hold close to the Brinkman-Rice transition on the 
metallic side. 

The proper form of the wave function for arbitrary band degeneracies (7) was first 
given by Biinemann and Weber || for translationally invariant systems [g^j = gi). 
Using the Gutzwiller approximation they derived the equations (39), ([48]), and (|49|) . In 
this work we showed that their results become exact in infinite dimensions. Furthermore, 
we explicitly covered all cases of symmetry breaking in the one-particle product wave 
functions |$o); the cases of spin and orbital ordering in materials with degenerate bands 
can equally be studied with the help of Gutzwiller-correlated wave functions. 

4-2. Slave Boson Mean-Field Theory 

For the one-band case the Slave Boson mean-field theory |J, fl3~l]-||33| yields the 



same results as Gutzwiller-correlated variational wave functions in infinite dimensions, 
including the cases of broken symmetry [15[ |16| . Therefore, the Slave Boson approach 



is variationally controlled at zero temperature in this limit. The Slave Boson saddle- 
point free energy can be re-derived from a "variational partition function" |15|, [T6|j . 
This construction shows that, at best, the Slave Boson mean-field theory is applicable 
up to excitation energies for which Fermi-liquid theory is valid. Since Fermi-liquid 



parameters can be derived from the variational ground-state energy [24|, the low-energy 



properties can equally well be described with the variational and the Slave Boson mean- 
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field approach. 

Recently, the Slave Boson mean-field theory for degenerate Hubbard models was 



worked out by Hasegawa [11] who used the extension of the Kotliar-Ruckenstein 



approach for the degenerate Anderson model by Dorin and Schlottmann ]pT0| . His results 
for N = 2 completely agree with ours and, thus, the above statements on the virtues and 
limitations of the Slave Boson mean-field approach also apply for the case of degenerate 



bands. Independently, Fresard and Kotliar [|T^] derived a set of Slave Boson mean-field 



equations and, as an application, reproduced Lu's results on the Mott transition 0. 

4-3. Anderson Transition in Strongly Correlated Electron Systems 

Our approach also covers the case of local, energetically random impurity potentials e? . 
For illustrative purposes we restrict ourselves to N = 1, U = oo, and no spin-flip 
hopping. In this case the system is a Mott insulator at half band-filling, 5 = 0, where 
6 = 1 — (1/L) J2i n i is the doping degree of the lower Hubbard band 0. In this case the 
variational ground-state energy becomes 

£ var ($o) = e c;y w^'n <I, o) + E ^ > ( 5l «) 

i^j;cr i 
%a = z • (516) 

1 - rti. a 

Recall that the local occupation densities n» = n^ a + n i; _ CT < 1 are given by n i]a . = 
(rii-a) = ($o|*vl$o)- 

We may interpret our variational results (51) in terms of a two-fluid picture 
which naturally arises for strongly correlated disordered electron systems. A recent 
introduction and overview on the theory of the Anderson transition in interacting 



electron systems is given in | 50fl . (i) A certain fraction of electrons may localize, e.g., 
a a electron is localized on site I in \$>o), n i — n i,a = 1- From the expectation value 
for the ground-state energy (51) it follows that the probability is zero that the electron 
hops off the site I. In addition, a —a electron will not hop onto this site because this 
is dynamically forbidden, qi- a = 0. Within our variational description the occupation 
of site I does not change. The site I is equally likely occupied by a a or a — cr electron 
and, therefore, the site contributes to a Curie-like susceptibility, (ii) The excluded 
sites represent an unretarded (random) hard-core potential for the remaining electrons 
besides the fluctuating local potentials and the dynamical constraint of no double 
occupancy. The strength of the impurity potentials and the degree of doping determine 
the fraction of "localized" and "mobile" electrons. 

For a small doping of the lower Hubbard band, 8 <C 1, even the "mobile" fraction 
of the electrons cannot carry DC current since their wave functions do not spread 
over macroscopic distances. The doping has to exceed some "percolation threshold" 
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to guarantee a finite DC conductivity. Above the critical doping for the Anderson 
transition 5 C one should observe a metallic conductivity at zero temperature but the 
magnetic signatures of the localized fraction of the electrons ( "local moments" ) should 
remain. It requires further doping beyond 5 C to destroy the Curie behavior of the 
magnetic susceptibility. 



Qualitatively similar results are found by Dobrosavljevic and Kotliar [34]. They 



used the dynamical mean-field theory approach and approximated the remaining local- 
impurity problem with the help of the Slave Boson mean-field approach of Barnes, 
Coleman, and Read and Newns [|3T|-|33|]. They also address the strong-coupling limit 
and find that local moments form in the metallic phase before the Anderson transition 
takes place. In their case the strength of the impurity potentials is varied for fixed 
doping. In both cases the position of the mobility edge is varied with respect to the 
Fermi energy, and the physics should qualitatively be the same. It should be clear, 
though, that our variational approach gives only a rather crude description for the local 
moment formation and the Anderson transition in the lower Hubbard band. 



5. Summary 

In this work we introduced a general class of Gutzwiller-correlated wave functions 
for multi-band Hubbard models. In contrast to earlier generalizations of the original 



Gutzwiller wave function |19[ to the case of degenerate bands [JIB], [2(], ^TJ we introduced 
independent variational parameters for each multiple occupancy of a lattice site. Only in 
this most general form each of the original variational parameters g^i can be expressed by 
the net occupancy densities my which give the average probability that configuration / 
is present on site i. 

We developed a diagrammatic formalism which allows the approximation-free 
evaluation of our general Gutzwiller-correlated wave functions in the limit of infinite 
spatial dimensions. Our analytical results reproduce recent results within the Gutzwiller 
approximation ||-|| and extend these to the case of a broken symmetry. In addition, 
the Slave Boson mean- field theory []10|— 112| becomes variationally controlled at zero 



temperature in the limit of infinite dimensions. As can be inferred from the single-band 
case, the Slave Boson approach at finite temperatures is limited to the Fermi-liquid 
regime [|15[ [l(|; see also Chap. 3.5]. Lastly, we briefly discussed the results of the 
variational approach to the Anderson transition in strongly correlated one-band systems. 

In principle, we could go beyond the limitations of the Gutzwiller approximation 
and systematically calculate ground-state correlation functions and 1/d corrections. 
However, the diagrammatic evaluation would be very tedious, and we consider it unlikely 
that such an effort will lead to qualitative changes of our understanding of band- 
degenerate Hubbard models based on our variational approach. For the moment it 
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appears to be more rewarding to elucidate the ground-state energy (|49|) in more detail. 
Thus far this expression was investigated only for rather simple model systems [frj- 



ljj . These preliminary studies showed significant differences for the ground-state 
magnetization as a function of the interaction strength between the Hartree-Fock 
and the Gutzwiller-correlated wave functions, as expected from the one-band case. 
Furthermore, discontinuous metal-insulator transitions occur for N > 2 for which the 



gap jumps to a finite value at the critical interaction strength || |TT . 

The investigation of the variational ground-state phase diagram for realistic models, 
e.g., for five d orbitals, is a numerically difficult task because the number of variational 
parameters exponentially increases with the number of orbitals. Nevertheless, one 
may use various symmetries between different multiple occupancies which considerably 
reduces the number of independent variational parameters such that the realistic case 
N = 5 for transition metals and their compounds should become tractable. 
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